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SUMMARY 

Dissemination prediction models have been reviewed to determine their 
applicability to a risk assessment for airborne carbon fibers. The review 
showed that the Gaussian prediction models using partial reflection at the 
ground agreed very closely with a more elaborate diffusion analysis developed 
for this study. For distances beyond 10000 m the Gaussian models predipted a 
slower fall-off in exposure levels than the diffusion models. This resulting 
level of conservatism was preferred for the carbon fiber risk assessment. The 
results also showed that the perfect vertical-mixing models developed herein 
agreed very closely with the diffusion analysis for all except the most stable 
atmospheric conditions. 


INTRODUCTION 

This study of dissemination models was made to support a risk assessment 
of carbon fibers released in aircraft accidents. 

Carbon fibers have light weight, high strength, and high stiffness and 
are being used together with epoxy resins to form composite materials. These 
materials may appear in large amounts in aircraft structures of the future. 

But the fibers have disadvantages, too. Their electrical conductivity has 
caused electrical equipment failures in the manufacturing process, and similar 
problems are likely to occur if the fibers are released from the composite. 

That could occur if the material is burnt in an aircraft accident. The risk 
study for that problem (ref. 1) involved the analysis of the atmospheric 
transport of the fibrous particles. The study reported here is a part of that 
effort; it represents a review of the existing dissemination models, their 
applicability, and some developments of simple models suitable for closed form 
mathematical analysis of such problems. 

The dissemination of lightweight particles and gases in the atmosphere has 
been studied both analytically and experimentally and involves disciplines 
such as fluid mechanics, meteorology, and statistics. An early summary text 
of the subject was prepared by Pasquill (ref. 2) . The main phenomena which 
must be accounted for in a dissemination analysis are turbulent mixing in the 
convection layer of the atmosphere, the downwind transport with the wind, and 
the falling of heavier-than-air particles. All of these processes are sto- 
chastic; therefore, the variables which must be predicted — such as concentra- 
tion at a point — have statistical distributions with wide variability. Only 
the long term averages can really be quantified by any simple mathematical 
model. Figure 1 shows a smoke plume from a fire test in which structural 
components made with carbon fibers were burnt so that the fibers could be 
released. Such a smoke plume can be idealized near its source as a conically 
expanding point source as shown by the solid lines. 



The simplest mathematical models that describe the particle dissemination 
are the Gaussian models such as the Turner model (ref. 3), and the Cramer models 
(ref. 4) . They are based on the fact that a point source of particles will 
produce such a conically growing particle cloud with a Gaussian transverse 
concentration profile. The growth is limited by the thin turbulent air layer 
between the ground and the inversion layer. The Gaussian models simply use a 
light beam reflection analogy to determine the growth of the Gaussian cloud 
beyond the point where the conically growing cloud first intersects either the 
ground or inversion planes. The falling of heavier-than-air particles and 
their deposition on the ground are simulated by making the ground reflectivity 
imperfect and by tilting the initial conical cloud at an appropriate fall angle. 

The more complex models are the particle-in-a-cell models (refs. 5 and 6), 
which solve the turbulent mixing, the downwind transport, and the falling in 
a three-dimensional analysis. The solutions can only be obtained numerically 
and require a large computer. 

The assessment of the effects of carbon fiber dissemination from aircraft 
fires was made by thousands of repetitive simulation of potential accident 
scenes for a number of cities around the nation. In these simulations the 
weather was randomly assumed to be characterized by one of the six Pasquill- 
Gifford stability classes (ref. 1) . The simple Gaussian models were computa- 
tionally ideal for this repetitive application. They have been developed to 
predict exposures close to the source of pollution and such predictions had 
been experimentally verified for distances up to 10 km. But in the carbon 
fiber study these models were to be applied up to distances of 100 km. One 
purpose of this study was to assure that the models were accurate, and that the 
empirical constants describing the particle deposition were correctly modeled 
for carbon fibers. 

This study recognized that there was little experimental data for diffusion 
at long distances and that because of the large variability of the atmospheric 
processes no single test series could be conducted to prove the accuracy of 
these models. Therefore, this study aimed to show analytically that these 
models are logically consistent and to determine if they satisfy the principle 
of mass conservations and an appropriate particle deposition law. 

For that step, a diffusion model was set up which handles the assumed 
deposition law and the diffusivity limitations at the ground in an exact manner. 
The predictions of the two Gaussian models was then compared with the diffusion 
model at points within a small distance from the source. Also, the predictions 
from all models had indicated that the particle distribution in the vertical 
direction become uniform between the ground and the inversion. Therefore, an 
analytical model was set up in which a uniform vertical particle distribution 
was assiomed. This model was used to show the parametric relation between the 
variables in the dissemination process and is useful for fast and simple expo- 
sure calculations. The derivation of all of the models are presented here and 
comparisons are made among the predicted results from the models. 
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diffusivity between adjacent cells 

integer coefficient 

total number of particles 

summation counter 

particle release rate 
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initial condition for particle flux, f/s 

constant in diffusivity equation 
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3 relative fall angle 

A finite increment for any function 

e thickness of an air cell, m 

Y vertical spread angle for rectangular cloud 

(j) (x) cloud expansion function 

0 declination of a point x below the source 

6^ azimuth spread angle 

9 elevation spread angle 

E 

T nondimensional function 

T initial condition for T 

o 

a dispersion in the crosswind direction 

y 

dispersion in the vertical direction 

0) width of a idealized cloud 

0) initial condition for cj 

o 


ANALYSIS 

Definition of Terms 

Because this study deals with particulate pollution, the quantities of 
pollutants are treated by enumeration, rather than by volume or mass. The 
particle concentrations are given in terms of particles per cubic meter and are 
given the symbol C. 

Many of the effects of atmospheric pollution are not just a function of 
the particulate concentration, but of the time-integral of concentration or the 
product of concentration and time. This quantity is termed either the expo- 
sure or dosage, has the symbol *E,* and is measured in particle seconds per 
cubic meter. 

Heavier-than-air particles may be tossed around by the turbulent atmo- 
sphere, but inside the turbulent cells of air the one constant force always 
acting on these particles is gravity, and it will move these relative to that 
air toward the ground at a mean fall velocity versus equal to that observed in 
still air. 
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Ultimately, all particles will deposit on the ground. The measure of 
ground contamination is the particle deposition density D measured in 
particles per square meter. The fact that all particles once airborne and 
contributing to the exposure must ultimately end up as a deposit leads to a 
mass conservation rule which all dissemination models should satisfy. 


The Mass Conservation Rule 

If particles fall with a mean fall velocity Vg, then they will deposit 
from a volume where the concentration is C at a rate CVg and, when all 
the particles have been deposited, the total deposition density will be 

D = jcv^ dt 


or 


D = Ev 

s 


If a total of N particles were initially airborne, the sum of all the depos- 
ited particles over all areas where particles had spread must be 

fo dA = N 


This is the mass conservation rule applicable to particle dissemination. If 
we combine this with the relation between exposure and deposition we get 

Je dA = N/v^ 

That relation prescribes a limit to the total area coverage which can result 
from a source of N particles having free fall velocity Vg. 

This rule should be satisfied by any dissemination model for heavier- than- 
air particles. 


The Gaussian Models 

For a continuous point source emitting particles at a rate n, the dif- 
fusion in the vertical and crosswind direction leads to a concentration profile 
with Gaussian distributions in these two directions. If z is the vertical 
distance and y the crosswind distance from the centerline, then the con- 
centration is 
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Because the functions are fully separable in y and z the concentration 
profile can always be obtained from the transverse or crosswind concentration 
integral (see appendix A) where 
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The reason for working with the crosswind concentration integral is that 
in most dissemination problems the dispersion in the crosswind direction is 
essentially unrestricted and the crosswind distribution can be assumed to 
remain Gaussian. In the vertical direction, the growth of the Gaussian plume 
becomes constrained by the lack of turbulent mixing into the inversion and the 
ground boundary layer. These complications are more easily solved by dealing 
with the crosswind integral I, which is independent of y. 


The accommodation of the zero turbulence both at the ground and at the 
inversion in the Gaussian models is approximate. Inplace of solving the 
diffusion equations, they assume that the particle behavior near the ground is 
analogous to the behavior of a light beam reflected off the ground. Figure 2 
shows a fiber source at a height h above the ground when the inversion is at 
a height H. The centerline of the "beam" is declined at an angle 3* If the 
reflectivity at the ground is designated as *r,* then at any point x downwind 
from the source and a distance Z above the ground, the contribution to the 
crosswind integral from the prime source is 


I 

x,z 


n 

o 

z 


-z^/2a 


2 

z 
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where 
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z = h - 3x - Z. 


Those Gaussian models which use imperfect reflection to model deposition 
(ref. 7) are referred to as the Gaussian partial reflection models. The reflec- 
tions of the prime source and the reflections of reflections contribute to the 
vertical distribution of crosswind integrals. It can be shown that if 

i = integer value of n/2 = INT (n/2) 
j = INT ((n - l)/2) 
k = INT ( (n + l)/2) 


uv^a 

Z 

where 

z^ = 2jH + h “ 3x + (-l)^Z 

and 

= 2kH - h + 3x + (-l)'^Z 


oo 

E 

n=l 


2,2 2,2 

. -zy2a . -zV2a 

1 1 z . X 2 z 

re + r e 


The above equations represent the Gaussian partial reflection model, in 
which the partial reflection coefficient r must be empirically determined 
to fit the particle deposition rate. The process is best understood by an 
illustration of a slowly growing light beam. 

Figure 3 is an isometric representation of the strength of the crosswind 
integrals at 10 downwind sections. The source beam has a Gaussian distribution 
defined by 


a 

z 


E 
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It is a linearly growing 'beam.' The reflection coefficient is unity both at 
the ground and at the inversion. The light beam alternately bounces off the 
ground and off the inversion. As it spreads the distortion from the reflections 
slowly changes the distribution from a Gaussian to a uniform distribution. The 
total section flux at all sections is 

F =u| Idz=ft 

J 

*'0 


as required for conservation of mass. In the particle dissemination problem 
the same physical reflection concept is applied to model the lack of turbulent 
mixing in the inversion and near the ground. The deposition of particles is 
modeled by partial reflection on the ground; and the fall velocity of the 
particles is modeled by inclining the cloud of an angle 3/ determined from the 
fall velocity and the wind velocity. 

For the small spreading angle shown in figure 3, this modeling would not 
very satisfactorily describe the behavior of particles because the reflection 
actually changes the direction of fall leading to a net vertical rise velocity 
for the particles remaining airborne after the reflection. This represents an 
anomaly in the modeling which shows up most strongly for particles with high 
fall velocities in very stable weather conditions. 

Unlike the example shown in figure 3, most particle dissemination problems 
have dispersion angles 0g larger than the fall angle 3- For those problems 
the vertical distribution soon becomes uniform and the bouncing light beam 
effect of figure 3 is not visible in the solution. Figure 4 shows a solution 
for a dissemination problem in stable weather with the following parameters: 


H = 100 m h = 50 m 

u = 1.5 m/s V =0.02 m/s 

s 

= 0.025 r = 0.3 

E 

At the last section shown, x = 2000 m, the distribution is almost uniform. 
The mean fall distance of the particles is 26.7 m as indicated by the solid 
centerline of the expanding plume. For the neutral and unstable weather catego- 
ries the dispersion angle 0g is even larger and the distributions become 
uniform very much faster. 

The location of the maximum exposures was important to the risk measurement 
and was obtained as an analytical expression here. It can be shown that the 
original source and its ground reflection contribute significantly to the 
maximum exposures at ground level. If the crosswind integral Ix,o consists 
only of these two terms, then 
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x,o 


= (1 + r) 
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2/ofl2 2 

-z /20gX 


/irr 6_x 
E 


in which 


z = h - 3x 


For this expression it can be shown that the maximum occurs at 


3 


x^ = ^ X 2G \Jl + G^-G 
Imax 3 J T 

and if the lateral dispersion of the cloud is linear the maximum concentration 
of exposure occurs at 

X = ~ V G I ^2 + G^ - g| 

The use of the partial reflection concept for simulating particle deposi- 
tion leads an anomaly. Conceptually we require that the deposition rate at 
the surface should be the product of concentration and fall velocity so that 
the crosswind integral of deposition rate should be The total deposi- 

tion between the source and section x is 
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•'0 


I V 
x,o s 


dx 
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and the residual flux is 



Mass conservation requires that the sum of these two integrals be equal to 
the particle flux n 


n = u| I_dZ+v I I 

I x,Z s I x,o 


dx 


The Gaussian models have a different relation governing deposition. This 
can. be seen at small distances x where only the original source and its 
first reflection contribute to the ground level concentration - 

Let 


I 

x,o 


(1 + r) 


n 

u/^ a 

z 


“ Z ^/20 


2 

z 


e 


where 


z == h - 


The flux through the surface is a vector with a direction shown in 
figure 5, where 


0 = tan 


1 


(h/x) 


The net vertical component from the flux of the prime source and its first 
reflection is 


F 
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(1 “ r) ~ u I 

X 


x,o 
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Because this deposition flux has a different functional form (x in 
the denominator) than the desired deposition rate VglxjO the Gaussian partial 
reflection model cannot satisfy the desired deposition rule, and the mass con- 
servation rule simultaneously. 

One alternative which had been used (ref. 8) is a Gaussian depletion 
model. That model uses full reflection of the particles, together with an 
overall depletion function such that 



and 

n -z^/2o^ 

I = — y 

u/^ 


plus reflection terms. 

The Gaussian depletion model satisfies the mass balance equation, and the 
desired deposition law. However, especially for heavy particles, the full 
reflection shows the full light beam characteristics and produces sometimes 
uneven, discontinuous particle distributions. 

The model described in reference 8 uses the Slade depletion functions. 

The comparisons in this report are based on the Gaussian exponential 
depletion model in which 


n 


n e 


-3x/h 


which is the depletion function obtained for the perfect vertical mixing models 
described in the following sections. 


A Differential Equation Diffusion Model 

The anomalies identified with the Gaussian dispersion models, and also 
the anomalies in the closed fom models are related to oversimplifications 
in the vertical direction. 

In this section a model is presented in which the applicable diffusion 
equation is first solved in the vertical direction between the ground and the 
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inversion. The solution in terms of the crosswind integral as a function of 
time can then be turned into a dissemination solution. 

The diffusion equation in which the diffusivity is nonconstant and in 
which the particles have gravitational settling is 


3t 


3^1 

K — + 
3z"^ 




9z 


A physical model resulting in this diffusion equation is presented in 
appendix A. The diffusivity can be a function of z and t. Because very few 
problems with variable diffusivities can be solved in closed form, the solution 
was formulated in a finite difference form. 


The determination of the diffusivity presented an anomaly. The solution 
of the differential equation for a point source in a field of uniform dif- 
fusivity is a Gaussian profile 
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uv^ a 
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-z^/20 
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where 


a = /2Kt 

Z 


The Gaussian models use this distribution profile, but consistent with 
observations generally stipulate that the vertical dispersion grows linearly 
with time or distance rather than with the square root of time 


a = 6„x or a 
z E z 


0_ut 

E 


To obtain solutions from the differential equation, which use the same 
rate of growth of the parameters 0^, the diffusivity was made a function of 
time. 

If we require a dispersion = ^E^^ ^ rate of growth 


da 
z 

dt 
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then 


a 

z 


da 

z 

dt 


0 u t 


da 

z 

In a previous equation, = K so that the diffusivity in a vertical 

plane at the downwind distance x will be 


K = Q^u^t or K = 6^ux 
r] £ 

This assumption will give the diffusion model the same characteristics as the 
Gaussian models as long as the cloud can grow without constraints. 

In modeling the convective layer the dif fusivities should go to zero both 
at the ground and at the inversion. For the purposes of this study the 
following function was chosen for the diffusivity. 


where p was arbitrarily chosen to be p = 0.1 for most calculations. 


The Perfect Vertical Mixing Model 

When the turbulent mixing occurs rapidly, and the fall velocity of the 
particles is small, their distribution in the vertical direction very rapidly 
becomes uniform. 

If perfect mixing in the vertical direction is assumed as an initial 
condition, simple closed form mathematical models can be derived which satisfy 
the requirement that the local deposition rate is the product of concentration 
and deposition velocity. Two idealized plume configurations have been treated, 
a vertical line source and a ground based point source. The vertical line 
source represents a good idealization for a fine plume which grows almost 
vertically and which gives off particles uniformly along its height. A ground- 
based point source represents a good idealization for a time plume in a high- 
wind in which the plume essentially remains attached to the ground- A simple 
model of a particle cloud from a continuously emitting source was created as 
shown in figure 6. The cloud has a rectangular cross section of height h and 
width 0) at a distance R downwind from the source. 

To conserve the mass in the cloud, the sum of the particle flux through 
section B and the particle deposition in the element area 0)Ax must equal the 
influx of particles through section A. 
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If the particle concentration C is uniform and the deposition velocity 
is Vg then the deposition rate in the element is 


C V U)AC 
s 

The outflux of particles is 


(C + Ac) (o) + Aco) (h + Ah)u 


and the influx of particles is 


C 0 )h u 


This results in a differential equation for the concentration profile 


dx \h h dx 0) dx/ 


where 


3 = is particle fall angle. 


If 


c() (x) 


h h dx 03 dx/ 


then the concentration profile is given by 


c = c 

If n is the fiber flux at z^, through a cloud h^, O)^ then 


C h OJ u = n 
o o o 

provides the input condition for the solution of the constant C* . 
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The vertical line source .- To model the vertical line source the height 
of a particle cloud is a constant H, and the width of the cloud is proportional 
to the radial distance. Therefore 


h = H and o) = ax 



The particle flux at x == x^ defines the initial condition so that 

C Hx u = n 
o o o 

which gives 

C = i 

Hau X 

The exposure profile is the time integral of concentration 

E = Jc dt 

and since 

N = Jn dt 

„ N 1 -3x/H 

E = e ^ ' 

Hau X 
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The local ground deposition density is 


D = E V 


so that 


n = JL_ e-3x/H 
Hax 


This satisfies the condition that the total deposition must equal the number 
of particles 


f 

•'n 


D dA = N 


It is important to realize that the crosswind integral of concentration 



is independent of the lateral spreading of the cloud. 

If, instead of assuming a uniform lateral concentration profile, we 
assume a Gaussian profile such that 
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whose crosswind integral is 


I = c /H o 

x,o y 


then 
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uh/^ a 
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SO that the concentration at any point is 
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x,y 


a uH 
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-3x/H -y^/2a^ 

e ^ 


and 
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x,y 


e e 

/2tt a uH 
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The ground-based point source To model a ground-based point source the 
height of the cloud is assumed zero at the source. If we assume the height to 
vary as 


h . h(i - 


and 


oj = ax 


We get a rectangular cloud growing from a point on the grour^d for which the 
expansion function 


({) (x) = 




-yot/H 


l(l - e-Y>‘/«) „(l - 


To simplify the expressions we substitute 


-Yx/H 

e ' = T 


The concentration profile is 


aHux \ / \1-T/\1-T/ 
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and the exposure profile is 



and the crosswind concentration integral 



If the crosswind concentration profile is Gaussian 


C 

x,Y 
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-y^/2a 


2 

Y 


where again the lateral dispersion can be chosen to fit the dissemination 

data. 


Most frequently/ predictions from the diffusion models are required 
in the form of contours of (constant) exposure. 

The equation for the distance y of the contour *E* from the centerline 
of calculations is 



Figure 6 shows a contour plot for the release conditions assumed for neutral 
atmospheres. 

N = 10^ H = 400 m 

= 0.02 m/s u = 6 m/s 

0^ = 0.07 
E 

= 0.07 
Y = 0.0877 
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MODEL COMPARISONS 


Elevated Sources 

Four models have been presented in the previous sections. They are the 
Gaussian partial reflection model, the Gaussian exponential depletion model, 
the perfect vertical-mixing models, and the analytical diffusion model. 

These four models will be compared in this section. All models have 
identical spreading in the transverse or crosswind direction because the 
almost linearly growing Gaussian dispersion represents the only existing 
data, and has been incorporated in all models. The comparisons are, therefore, 
made in terms of crosswind integrals. 

All four models were used to calculate the crosswind integrals for ten 
downwind sections and all points between the ground and the inversion. 

These calculations were made for two particle fall velocities, 0.02 m/s 
and 0.2 m/s in three meteorological conditions. 

Table I shows a matrix of model parameters for three meteorological 
conditions . 


TABLE I.- MODEL PARAMETERS 


Condition 

Inversion 
height, H 

wind 
vel., u 

Vert. dlsp. 
angle, 0^ 

stable 

100 

1.5 

0.025 

Neutral 

400 

6.0 

0.07 

Unstable 

1000 

2.0 

0.15 


The results for the vertical particle distribution are plotted as isometric 
views of the crosswind integrals at ten downwind stations in figures 8 to 13. 

Table II shows the matrix of conditions applicable to each of these 
figures. 
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TABLE II.- MODEL INPUT CONDITIONS 


Meteorological 

category 

Particle velocity 
m/s 

Isometric 

figure 

Ground- level 
figure 

stable 

0.02 

8 

14 

Stable 

0.2 

9 

15 

Neutral 

0.02 

10 

16 

Neutral 

0.2 

11 

17 

Unstable 

0.02 

12 

18 

Unstable 

0.2 

13 

19 


While figures 8 through 13 show the distribution of the particles in the 
vertical and downwind directions, a more precise comparison of the models is 
presented in figures 14 through 19. They are semi logarithmic representations 
of the ground-level values of the crosswind integrals. 

The isometric sections show that in five of the six conditions, the 
turbulence is strong enough to produce essentially perfect vertical mixing. 

The isometric sections show that in neutral and unstable weather 
(figs. 9 through 12) the vertical distribution becomes uniform before signifi- 
cant depletion occurs, even for the high fall velocity. In the stable atmo- 
sphere and the 2 cm/s fall velocity the distribution becomes triangular 
(fig. 8) at the 2000 m range, and for the 20 cm/s fall velocity (fig. 9) the 
cloud falls to the ground before dispersing significantly. 

The ratio of dispersion angle Gg to the fall angle 3 is the logical 
measure of the dispersive strength. From these six cases we obtain the result 
that the cloud will fall to the ground if Gg/B < 1, and that it will deplete 
after distributing uniformly if Gg/3 >2. 

The semilogarithmic presentation of the ground-level values of the cross- 
wind integrals show the smaller differences between the four models in detail. 
Figure 14 shows the solutions for the 2 cm/s fall velocity particles in the 
stable atmosphere. The straight solid line is the result for the perfect 
vertical-mixing model. The Gaussian exponential depletion model rapidly 
approaches this solution. The analytical diffusion model also assymptotically 
approaches this solution. Only the Gaussian partial reflection model which was 
tuned to agree with the analytical diffusion model 2000 m downwind from the 
source deviates from the other model results. This trend agrees with the 
analytical observation that the mathematical form of the deposition by partial 
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reflection does not agree with the constant fall velocity deposition. The 
same trend is visible for the other atmospheric conditions, especially in 
figures 15, 16 and 19. In all cases, the reflection coefficients had been 
determined by iteration so that the Gaussian partial reflection model would 
agree with the analytical diffusion model at or just beyond the peak values 
of the crosswind integral. Table III lists the values of the partial reflec- 
tion coefficients used in the six analyses. 


TABLE III.- PARTIAL REFLECTION COEFFICIENTS 


Condition 

Light part. 
2 cm/s 

Heavy partic. 
20 cm/s 

Stable 

r = 0.5 

r = 0,01 

Neutral 

0.85 

0-45 

Unstable 

0.9 

0.35 


While the use of the partial reflection coefficient allows tuning the model 
for a good agreement at one distance, it is important to realize that, at 
long distances downwind, the model will predict conservative values of con- 
tamination . 

The Gaussian exponential depletion model in all cases predicted the 
highest peak values and converged to the solution of the perfect vertical 
mixing model. In the stable atmosphere with the heavy particles this solution 
(fig, 15) breaks down. It converges to the solution for the perfect vertical 
mixing model, although because of the low turbulence, this weather situation 
never produces perfect vertical mixing for such heavy particles (fig. 9) . 

Also, the bouncing light beam effect is very visible in this solution. Not- 
withstanding those shortcomings, the solutions predict the correct maximum 
values, and the errors only become large when the exposures are very much 
lower than the maxima observed. 


Ground-Based Point Sources 

The perfect vertical mixing model with the ground-based point source has 
been compared to the solutions obtained from the diffusion model. The eleva- 
tion growth angle of the cloud y was assiamed to be equal to the elevation 

dispersion angle 0^ x A/2 to match the maxima of the original Gaussian and 
uniform distributions - 
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Figure 20 shows the semi logarithmic presentation of the ground-level 
values of the crosswind integral. The two models agree closely over the 10 km 
range. Similarly, good agreement was obtained for neutral atmospheres 
(fig. 21) and the unstable atmosphere (fig. 22) . The ground-based point source 
model has the great advantage that contamination levels can be computed from 
a single closed-form algebraic expression, while the Gaussian models contain 
series of terms from successive reflections, and the diffusion model requires 
a slow numerical solution of its differential equation. 


CONCLUDING REMARKS 

The Gaussian dissemination models which were used in the carbon fiber 
risk assessment, have been compared to an analytical diffusion model. Partial 
reflection coefficients were obtained for pollutant particles of 2 cm/s and 
20 cm/s fall velocities, characteristic for single carbon fibers and small 
clumps. The model predictions based on these reflection coefficients were in 
close agreement with the solutions from the diffusion model for similar atmo- 
spheric conditions. For large downwind distances the Gaussian models tended 
to predict higher than expected exposures. For the graphite fiber risk assess- 
ment, the use of these Gaussian models is expected to provide reliable, if 
somewhat conservative, exposure calculations. 

Because most pollutant cloud models showed that uniform vertical mixing 
over the convective air layers was reached very rapidly, a perfect vertical 
mixing model was established which is useful for fast and simple exposure 
calculations. Those models provided a parametric analysis technique for the 
graphite risk assessment. 
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APPENDIX A 


DIMENSIONAL CONDENSATION OF THE DIFFUSION PROBLEM 

For a continuous particle source, the concentration of particles in a 
cell drifting with the wind can be expressed by a two-dimensional diffusion 
equation 


at y 8z ^8z ^sj ^z 


( 1 ) 


There are some solution types for this equation which have the form 


C = C F^(y,t)F2(z,t) (2) 

in which the functions of y and z are separable. If the function describ- 
ing the distribution of the particles in the crosswind direction, Fj^(y,t), 
has the characteristics of a distribution function, such that its integral in 
the y direction is a constant. 



F^(y,t) 


dy = B 


( 3 ) 


then the two-dimensional diffusion problem can be condensed into a one-dimen- 
sional problem in the vertical direction. 

The various differentials of the concentration profile are 


3^- C'Fi(y.t) ^Pi(y,t) 

C-F^Cy,t) ^F^(z,t) 
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.2 .2 

— . C'Fj(y,t) — F^Cz.t) 

dz dz 


2 2 

CF^Cz.t) ^Pi(y.t) 

dy dy 


When these are substituted into equation (1) 


F, (y,t) F^{z,t) + F^(z,t) F^ (y,t) = KF^ (y,t) F^(z,t) 


9t "1’ 


9z 


2 2 


+ F^(y,t) -f F^lz.t) 


+ KF2(z,t) F^{y,t) 

3y 


(4) 


If this equation is integrated with respect to y 


■^F^Cy^t) I F^(y,t) dy + F^Cz^t) j ^F^(y,t) dy = K — - F^(z,t) 

CX> •'—00 ^ ^ 


>= f F^(y,t) ay + 

oo 


f 

•'—O 


f iIf 


X 1 F^(y,t) dy + KF 2 (z,t) 


(y,t) dy 


(5) 
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Since F^(y,t) is a distribution function which has 


OO 


t) = B 


then 


( F^(y,t) dy = f 

cx> •/ o 


• 00 2 

(y,t) dy = I F^(y,t) dy = 0 

3y 


When these are substituted in equation (5) , we get 

B ^ = BK ^ Fjlz.t) + b(|S + A 


( 6 ) 


If I is defined as the crosswind integral of concentration 

(y,t) dy 


I = C'F^{z,t) 


^CO 

•'-.oo 


Then, equation (6) becomes 


3 3 

^ I = K — I + 
9z^ 


/3k \ 9i 

y9z '^sy 9z 


(7) 


This is a one-dimensional diffusion problem in the vertical direction. The 
dissemination models derived herein are first obtained in terms of the cross- 
wind integral. Any lateral distribution such as the Gaussian can be given to 
the crosswise concentration profile. 


26 



APPENDIX B 


A PHYSICAL MODEL FOR THE DIFFUSION EQUATION 

The atmosphere can be subdivided into a series of layers. Perfect 
mixing can be assumed within each layer together with a fixed air exchange 
between adjacent layers. 


Figure 23 shows a series of five layers of depth £ and lateral dimensions 
Ax and Ay. The layers are joined by imaginary airflow ducts exchanging air 
at a flow rates q. Simultaneously, the particles can fall through the layers 
at a fall velocity v^. 

Let the concentration around the center element be described by the 
Taylor expansion 


-o ^ If - i s 

dZ 


and let the flow rates q also be expressed as 


q(e) = q + £ + %- — § 


^o dz 


2 


If the concentration is defined as C(e) then 

c(e) =c =c + e|^ + ^ ^ 

+1 o dz 2 . 2 

dz 


and if the exchange flow between C and C is defined as 

o +1 


q(e/2) = q =q + ^ ^ 

' ^ 0,1 ^o 2 9z 8 .2 

dz 


The particle flow AN into the center element is then 


an = - Cq) + q^^_^ 
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The settlement flow through the boundary is 

At|v^ Ax Ay(c^j - 

After expanding concentrations and flow rates 

A., A.-) 2 a^c . 2 aq ac a a 9 

an = Atjq^e — + e ^ ^ Ax Ay e ^ 


ac' 

z I 


if we define an exchange flow velocity such that 


V Ax Ay = q 
a ^ 


then 


An ^ 3^ ^ ^ ^ 

At Ax Ay "^a ^^2 9z 9z ^s 9z 


If the product of layer height and exchange flow velocity is K = v e, 
then this equation reduces to the normal diffusion equation 

^ ^ 9^C . ^ ^ ^ 

9t 2 9z 9z ^s 9z 

dz 


The finite element formulation for this equation is 
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For the outer layers the equation will be at the upper bound 



and for the lower bound 



Those formulations provide mass-conservative solutions to the one-dimensional 
diffusion problem. 
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Figure 5. Definition of the Flux Vector 



Figure 6. Configuration of a Typical Section Through 
an Expanding Rectangular Cloud. 
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Figure 7. Exposure Contours for Ground-Based Point Source 
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in Neutral Atmosphere for Release of 10 Particles, 
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Figure 14. Variation of Crosswind Integral with 

Downwind Distance in Stable Atmosphere 
2 cm/s fall velocity 
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Figure 15. Variation of Crosswind Integral 

with Downwind Distance in Stable Atmosphere 
Fall Velocity 20 cm/s 
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Figure 16. Variation of Crosswind Integral 

with Downwind Distance in Neutral Atmosphere 
Fall Velocity 2 cm/s 
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Figure 17, Variation of Crosswind Integral 

with Downwind Distance in Neutral Atmosphere 
Fall Velocity 20 cm/s 
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Figure 18. Variation of Crosswind Integral with Downwind 

Distance in Unstable Atmosphere, 2cm/s Fall 
Velocity. 




Figure 19. Variation of Crosswind Integral with Downwind 
Distance in Unstable Atmosphere 20 cm/s Fall 
Velocity . 
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Figure 20. 
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Unit Crosswind Integral for Ground-Based Point 
Source in Unstable Atmosphere and 2cm/s Fall Velocity. 
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Figure 21. Unit Crosswind Integral for Ground-Based Point 

Source in Neutral Atmosphere and 2cm/s Fall Velocity 



Figure 22. Unit Crosswind Integral for Ground-Based Point 

Source in Stable Atmosphere and 2cm/s Fall Velocity 





Figure 23. Finite Element Model for Atmospheric Diffusion 
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